home *** CD-ROM | disk | FTP | other *** search
/ QRZ! Ham Radio 8 / QRZ Ham Radio Callsign Database - Volume 8.iso / pc / files / dsp / 96000tar.z / 96000tar / 96000 / appb / b139.asm < prev    next >
Assembly Source File  |  1992-04-28  |  1KB  |  35 lines
  1. ; This program was originally published in the Motorola DSP96002 Users Manual
  2. ; and is provided under a DISCLAIMER OF WARRANTY available from Motorola DSP
  3. ; Operation, 6501 William Cannon Drive West, Austin, Texas 78735-8598.  For
  4. ; more information, refer to the DSP96002 Users Manual, Appendix B, DSP
  5. ; Benchmarks.
  6. ;
  7. ; B.1.39    Nth Order Polynomial Evaluation for Two Points 
  8.  
  9. ;An Nth order polynomial c1XN + c2XN-1 + ...cNX + cN+1 can be factored 
  10. ;and represented as ((c1X + c2)X + c3)X + ...) + cN+1. This routine 
  11. ;evaluates the polynomial at both X = s and X = t. 
  12. ; 
  13. ;Memory Map :   X             Y 
  14. ; 
  15. ;       r1 ->   s             t 
  16. ;               . 
  17. ;               . 
  18. ;       r0 ->   c1 
  19. ;               c2 
  20. ;               c3 
  21. ;               . 
  22. ;               . 
  23. ;               cN+1 
  24. ;  Setup N equ  order of polynomial 
  25. move #coef,r0 
  26. move #2_pts,r1 
  27. move x:(r1)+,d5.s y:,d4.s                     ; s, t 
  28. move x:(r0)+,d1.s                             ; c1 
  29. move d1.s,d0.s 
  30. ; Inner loop for evaluating 2 consecutive points 
  31. do #N,_loop 
  32. fmpy.x d1,d5,d1  x:(r0)+,d2.s                 ; c(n)*s, c(n+1) 
  33. fmpy   d0,d4,d0  fadd.x d2,d1                 ; c(n)*t, c(n)*s+c(n+1) 
  34. fadd.x d2,d0                                  ; c(n)*t+c(n+1) _loop 
  35.